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1 Abstract: 

A quantum mechanical approach is proposed for the for- 
mation of antihydrogen (H) in the ground and excited states (2s, 
2p) via the mechanism of three body recombination ( TBR ) in- 
side a trapped plasma of anti proton ( p ) and positrons ( e + ) 
or in the collision between the two beams of them. Variations of 
the differential ( DCS ) as well as the total ( TCS ) formation 
cross sections are studied as a function of the incident energies of 
both the active and the spectator e + s. Significantly large cross 
sections are found at very low incident energies in the TBR pro- 
cess as compared to other processes leading to antihydrogen. The 
present ( H ) formation cross section decreases with increasing 
positron energy ( temperature ) but no simple power law could 
be predicted for it covering the entire energy range, corroborat- 
ing the experimental findings qualitatively. The formation cross 
sections are found to be much higher for unequal energies of the 
two e + s than for equal energies, as expected physically. 

2 Introduction: 

Production of antihydrogen, the simplest and the most 
stable bound state of antimatter is one of the current topics nowa- 
days both from the experimental and the theoretical perspectives 
mainly because its study provides various fundamental differences 
between matter and antimatter. Particularly, cold antihydrogen 
(H ) atom is an ideal system for studying the fundamental sym- 
metries in physics e.g., the CPT invariance theorem in the stan- 
dard quantum field theory and the gravitational weak equivalence 
principle for antimatter. The major challenge facing the H re- 
search is the production of cold and trapped ground state H that 
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is needed for the precise laser spectroscopy. Apart from these, 
there are many important practical applications of the H out 
of which the followings are worthy to be mentioned. First, the 
antihydrogen may also be used for igniting inertial confinement 
fusion pellets, the feasibility of which was already investigated 
[1]. Second, the antihydrogen finds important applications in the 
propulsion system [ 2 ]. 

In view of the recent technological advances in the 
cooling and trapping mechanism of antiprotons ( p ) and positrons 
( e + ), the long term goal for the production of cold and trapped 
H, necessary for the high precision spectroscopic studies has now 
become possible. This has motivated theoretical workers to ven- 
ture different processes producing antihydrogen. The most im- 
portant of these processes is the following three body recombi- 
nation ( TBR ) p + e + + e + — > H + e + ( I ) in which a 
spectator particle carries away the excess energy and the mo- 
mentum released in the recombination. The above reaction poses 
to be more efficient by orders of magnitude [3,4] compared 
to other H production processes, e.g., the radiative recombina- 
tion ( RR ) [ 5 - 7 ] , the three body charge transfer between the 
Ps and the p [ 8 - 19 ]. The main reason for this is due to the 
following. The spectator positron in the TBR process efficiently 
carries off the extra energy, unlike the other H production reac- 
tions. Another important advantage of the process ( I ) is that 
the reactants are stable charged particles which can be held in a 
trap for cooling and then subsequently for the recombination to 
occur. In fact, it is found experimentally [3,4] that the TBR 
in the trapped plasma of antiproton and positrons happens to be 
the most efficient H production reaction at low and intermedi- 
ate energies. However, the main disadvantage of the TBR is that 
the H is favourably formed in the excited states [ 3 ], although 
for the high precision spectroscopic studies, the ground state H 
is highly needed. In contrast, in the RR process, although the 
ground state is favoured, the cross-section itself is much lower 
[3,4] . 



2 



Regarding the experimental situation for the TBR pro- 
cess, the three main International Groups are working on it with 
much endeavour at Cern, e.g., ATHENA [ 20 - 25 ], ATRAP 
[ 26 - 30 ] and Harvard [ 31 - 33 ] while another group from Riken 
[ 34 , 35 ] is also concentrating on the experiments of cold and 
trapped H production. In all the experiments attention is being 
paid mainly to the temperature dependence of the H production 
at extreme low energies in the range of meV. 

As for the theoretical situation, the first detailed study 
for the TBR is due to Robicheaux [ 36 ] in the framework of 
classical trajectory Monte Carlo ( CTMC ) method. However, in 
this calculation the Author introduced some fraction of electrons 
along with the e + as well as a strong magnetic field in order 
to make the process feasible. It was noted [ 36 ] that the H 
formation reduces substantially in presence of the e~s. Prior to 
and also following this work [36] , there exist some calculations by 
the same Author [ 37 , 38 ] that mainly study the temperature 
dependence of the TBR process based on some statistical models. 

The present work addresses the study of the H for- 
mation cross sections ( both differential and total ) through the 
TBR mechanism in the collision between the positron and the p 
plasma. To our knowledge, this is the first quantum mechanical 
attempt along this direction in the TBR process. Although ex- 
perimentally the TBR process favours the H formation in highly 
excited states, for the theoreticians it is much easier to calculate 
the cross sections in the ground and low lying states. Thus in the 
absence of any experimental cross section data, the present the- 
oretical estimates of the ground and excited ( 2s , 2p ) states H 
formation cross-sections might give some stimulus and guidelines 
to the future detailed experiments. 

The present model corresponds to the following situ- 
ations . We consider an ensemble of weakly correlated positrons 
and the e + plasma density is assumed to be low enough so that the 
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e + — e + interaction can be treated as a perturbation. The antipro- 
ton is treated as a stationary ionic target located at the origin 
of coordinates which corresponds to the experimental situation 
of a cold and trapped p . The latter assumption should be quite 
legitimate when the positron velocity is much faster than that of 
the antiproton which happens to be the case due to large mass 
difference between the two. Since the recombination reaction re- 
quires a third body for the energy and momentum conservation 
of the process, another e + ( electron ) of the plasma serves this 
purpose and the process becomes a TBR one. Further, there is 
a probability of exchange between the active and the spectator 
positrons which is also incorporated in the present model. 



3 Theory: 

The present study deals with the following three body 
recombination process : 

p + e+ + e+ -> H + e+. ( 1 ) 

In the present formulation the p is assumed to be stationary 
( target ). The prior form of the transition amplitude for this 
process is given by 

T if = (*7(fi, f 2 )(l + P) | Vt I Vi (ri, f 2 )> , ( 2 ) 

where P denotes the exchange operator corresponding to the in- 
terchange of the positrons in the final channel. Vi in equation 
( 2 ) is the initial channel perturbation which is the part of the 
total interaction not diagonalised in the initial state and ipi is 
the corresponding asymptotic wave function . The final channel 
wave function satisfies the three body Schrodinger equation 
obeying the incoming wave boundary condition; 

(H - E)V] = 0. ( 3 ) 

The total Hamiltonian ( H ) of the system can be written as 

ff = _!v5-i V i-i-i+i (4) 
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where f\ and represent the position vectors of the active e + ( to 
be transferred ) and the spectator e + s respectively. The atomic 
unit ( a.u. ) is used throughout the work. 

The initial channel asymptotic wave function ^ in equa- 
tion ( 2 ) satisfies the following Schrodinger equation 

(_i V ? -|Vl -i -\-E ) ^ = (5) 

and is given by 

^. = Nj e ik ^ 1F1 [ iaij, 1, -i {kfj - kj ■ r})], (6) 

with Nj = exp(— T(l—iaj);] = 1, 2 ; otj = — ^- ; kj denotes the 
incident momentum of the active or the spectator e + respectively. 
The approximated final state wave function is chosen in the 
framework of eikonal approximation as follows : 

*7 = Mn) J*'* exp [irjf (i 12 - \ 2 )dz'] (7) 

with rjf = kf being the final momentum of the spectator e + ; 

represents the bound state wave function of the H atom . 
Finally the differential cross section for the process 
( 1 ) is given by 

& = & [l(l/ + ff| 2 ) + f(l/-fll 2 )], (8) 

where f and g corresponds to the direct and the exchange ampli- 
tudes respectively. 

Using the following contour integral representations of 
the eikonal phase factors [ 39 ] as well as the coulomb functions 
[ 40 ] and after much analytical reductions [ 41, 42 ] , the transition 
matrix element ( 2 ) is finally reduced [ 43 ] to a three dimen- 
sional integral which is evaluated numerically by using different 
quadrature methods . The eikonal phase factor is of the form 

V± {i *- n) = snswfeiFWS) fri-^"- 1 exp(-\y)dy ( 9a ) 

where the contour c has a branch cut from to oo [ 39 ] ; the 
confluent hypergeometric function: 
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1F1 (ia, 1, z) = ^ jj° > p (a, t) exp(zt)dt ( 9b ) 

with p (a,t) = t~ l+l0i {t — l)~ m , Ti is a closed contour encircling 
the two points and 1 once anti - clockwise [40]. At the point 
where the contour crosses the real axis to the right side of 1, arg 
t and arg ( t-1 ) are both zero. 

4 Results and Discussions: 

We have computed the H formation cross sections both 
differential and total for the TBR process ( 1 ) in the frame- 
work of the coulomb distorted eikonal approximation ( CDEA ), 
where distortions have been included in both the channels. The 
exchange between the active and the spectator e + s is also incor- 
porated. Since the present process ( 1 ) is an exothermic reaction 
it can occur even at zero incident energy. However, our results are 
not converged below 5 eV due to computational problems and are 
therefore not reported here. Furthermore, it may be mentioned 
that the present model might not yield very reliable results at 
extreme low energies. 

Figs ( 1 - 5 ) exhibit the present differential cross sec- 
tions ( DCS ) in the ground and excited states ( 2s and 2p ) for 
different incident energies of both the positrons e.g., E\= E^ = 
10, 20, 25, 30 and 50 eV respectively. The figures reveal that 
the H formation ( in all the states ) is strongly favoured in the 
forward directions and as such the DCS are presented upto 60° 
only, beyond which the cross sections become negligible. At very 
low incident energy, the magnitude of the formation cross section 
is found to be largest for the 2p state and smallest for the Is state 
while the 2s lies in between, i.e., 2p > 2s > Is ( vide fig. 1 ). 
This trend of the DCS is noted upto 15 eV ( not shown in figure), 
although with increasing incident energy, the maximum of the 2p 
DCS decreases and tends towards the 2s maximum so that at E\ 
= E 2 = 20 eV, the 2s overtakes the 2p ( vide fig. 2 ). The DCS 
peak in this case is in the order 2s > 2p > Is. In contrast, at 
intermediate and high incident energies ( ~ 25 eV onwards ), the 
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DCS is maximum for the Is state and minimum for the 2p state 
while the 2s lies in between ( i.e. Is > 2s > 2p, figs. 3 - 5 ). 

As for the position of the DCS maximum, at low inci- 
dent energies ( ~ upto 15 eV ), the Is and 2p maxima lie at some 
lower scattering angles ( ~ 20° ) while the 2s maximum occurs 
at the extreme forward ( ~ 0° ). With increasing energy, the 
DCS maxima for these two states ( Is and 2p ) move towards the 
extreme forward ( ~ 0° ), while the 2s maximum moves in the 
reverse direction ( vide figs. 2 - 4 ). However, at high incident en- 
ergies (e. g., ~ 40 eV onwards ), all the partial DCS maxima are 
finally peaked at extreme forward 0° ( vide fig. 5 ), as expected. 

Figs. 6(a) - 6(c) again exhibit the partial DCS but for 
some unequal energies of the two incident e + s ( i.e. E\ ^ E2 ) 
along with a case for E\ = E2 ( 15 eV ) for the sake of comparison. 
The following interesting features are noted from the figures. All 
the partial DCS are found to be much higher ( by a factor of ~ 
2 to 2.5 ) when the energy of the active e + ( E\ ) is greater than 
that of the spectator one ( E2 ), i. e., when E\ > E^. The DCS 
for unequal ( E\ 7^ E2 ) energies lie much above than those for 
equal energies ( E\ = E2 ) . This is quite expected physically due 
to strong repulsion between the two e + s at equal energies. 

Next we come over to the total cross sections ( TCS ) 
for the H formation displayed in figs. 7 and 8 for different sets 
of incident energies. Figure 7 displays the partial TCS when the 
two incident e + s share equal energy (E\ = E2) . As in the case of 
DCS, at low and intermediate incident energies (~ upto 20 eV), 
the partial TCS follows the order 2p > 2s > Is ( inset of fig. 7 ), 
while beyond 25 eV it is in the decreasing order with excitation 
of the H state, i.e., Is > 2s > 2p. The dominance of the 2p 
TCS at low incident energies could probably be attributed to the 
long range polarization effects which is much stronger for the 2p 
state than for any other states. In fact, a major contribution to 
the polarization effect that mainly dominates at lower incident 
energies, comes from the lowest lying p state ( i.e., 2p state ). 
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Fig. 7 also indicates that although all the partial TCS 
decrease monotonically with increasing incident energy, they do 
not follow any simple power law ( e.g. ~ T~ 9 / 2 ), corroborating 
the experimental findings [ 21 , 23 ] . 

Figures 8(a) - 8(c) display the partial TCS against the 
active e + energy ( E\ ) for some fixed values of E% ( spectator ) 
while the insets exhibit the reverse one, i.e. TCS vs E<i for fixed 
E\. As in the case of DCS, for a fixed sum of E\ and E2, the 
partial TCS is found to be larger when E\ > E2 than for E\ < E^. 
Further, the TCS against E2 falls off much more rapidly than the 
TCS versus E\ ( cf. figs. 8 with their insets ). Regarding the 
relative magnitude of the partial TCS, for lower energy of the 
spectator e + , e. g., E2 = 10 eV ( fig. 8(a)), the general trend of 
the TCS follows the order 2p > 2s > Is at low and intermediate 
Ei while at higher E\, the above order changes to 2s > 2p > Is. 
Similar behaviour is noted in fig. 8(b) ( for E2 = 15 eV ) as in 
fig. 8 (a) with some exceptions at higher E\ . At intermediate 
E2 (25 eV, fig. 8 (c)), the partial TCS follow different orders for 
different ranges of E\, e.g., at lower E\ the 2p dominates while 
at higher E\ the Is dominates. However, at higher E2 ( ~ E2 > 
50 eV ), the Is cross section dominates through out the range of 
Ei except at very low energies ( E\ ~ 5 — lOeV ) where the 2s is 
most prominent ( not shown in figure ). 

Figure 9 demonstrates a comparison between the present 
TCS and the corresponding experimental results of Merrison et al 
[ 44 ] for the Ps(ls) +P — ► H(ls) + e + process. The experimen- 
tal Ps energies are converted to the antiproton energies following 
the relation : Ep(Kev) = ^^q* 918 . As was anticipated [3,4], 
the present TBR cross sections are found to lie much above the 
experimental data [44] for the abovesaid process. 

For the sake of some numerical measures, we have dis- 
played in Table I the present partial ( Is, 2s, 2p ) TCS along with 
some other existing theoretical results due to Mitroy et al 
[ 14, 17 ] for the process p + Ps — ► H(n, /, m) + e + using unita- 
rized Born approximation [ 17 ] and close coupling approximation 
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[ 14 ]. Results due to Sinha et al [ 19 ] for the above process in 
the eikonal approximation both with and without laser field are 
also included in the Table I. The incident energies are chosen in 
accordance with their [ 14 , 17 ] calculations. The field assisted 
( FA ) results [ 19 ] are presented for the field strength 0.01 a.u. 
and the frequency 0.043 a.u. 



Table I 



Energy 
(eV) 


Mitroy et al [14 ] 
Is ( 2s+2p) 


Mitroy et al [17 ] 
Is 2s 2p 


Sinha et al [ 19 ] 
FF ( Is ) FA( Is ) 


Present results 
Is 2s 2p 


13.60 


1.923 8.44 


1.46 0.66 3.81 


0.92 8.32 


373.07 705.4 530.06 


20.40 




0.94 0.254 1.76 


0.47 3.99 


68.25 59.62 36.13 


25.84 


0.735 1.729 




0.29 2.05 


23.16 15.43 11.048 


34.00 




0.394 0.08 0.396 


0.16 0.65 


5.14 3.29 1.21 


43.52 


0.24 0.28 




0.08 0.19 


1.38 0.73 0.21 


54.40 




0.13 0.03 0.06 


0.04 0.09 


0.507 0.212 0.651 


63.92 


0.078 0.053 




0.02 0.05 


0.165 0.064 0.034 



Table I again confirms ( as in fig. 9 ) that the present 
TBR cross sections are much larger than all the other processes 
leading to antihydrogen throughout the energy range considered. 
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Table II displays the probable power laws obeyed by 
the partial as well as the sum TCS for different incident energy 
ranges of the e + corresponding to figures 7 ( E\ = E 2 ) and 8 
(E 17 ^E 2 ). 



Table II 



Energy Range 

( in eV ) 


Power law obeyed 


7 — ' 7 — ' 

hi\ = h/2 


Is 


2s 


2p 


Is + 2s + 2p 


5 - 10 

10- 25 
25 - 50 


~ E~ 5 


~ E~ 37 
~ E- 5 - 3 




~ E~ 5A 
~ £T 5 - 9 


E 1 <E 2 { 10 eV ) 




5 - 10 


~ Er 1 - 6 


- Er 1 - 3 


-Er 1 - 8 


- Er 1 - 6 


E 1 >E 2 ( 10 eV ) 




10- 25 
25 - 50 


~ Ei~ l - b 

~ Er 1 - 5 


-Er 1 - 7 
-Er 1 - 8 


-Er 2 - 1 
-Er 2 - 3 


-Er 1 - 9 
-Er 1 - 9 



As is revealed from the table , the low energy partial 
TCS ( e.g. , Ei = E 2 ~ 5 - 10 eV ) falls off much slowly as 
compared to the intermediate and high energies and the slope of 
the Is TCS ( see also fig. 7 ) is much less than that of the others 
( 2s, 2p and Is + 2s + 2p ). Another important feature should 
be noted from this Table that for E\ ^ E 2 , the power of the 
exponent decreases as compared to the E\ = E 2 case throughout 
the energy range. This again indicates the better efficiency of the 
H production for unequal energies ( E\ ^ E 2 ) of the active and 
the passive e + s over a wider energy ranges ( E\). 
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5 Conclusions: 



The salient features of the present study are as follows: 

At very low incident energies the present TBR cross 
section for the H formation in the 2p state is found to be the 
dominant process among the three states Is, 2s, 2p while at inter- 
mediate and high incident energies, the ground state ( Is ) cross 
section dominates for both equal and unequal energies of the two 
positrons with some exceptions for the latter case ( E\ > E 2 ) . 

Substantially high cross sections are noted in the TBR 
model than in the other RR / charge transfer processes leading 
to antihydrogen. 

The partial TCS is found to be significantly higher 
when the active e + energy is greater than that of the spectator 
e + ( i. e., E\ > E 2 ) than for E\ = E 2 or for E\ < E 2 . 

For a more efficient production of H for a wider energy 
range, the unequal ( Ei > E 2 ) distribution of energy between 
the active and the spectator positrons could be suggested rather 
than the equal one ( E\ = E 2 ) . 

The present H formation cross section decreases with 
increasing e + energy ( i.e., temperature ) but does not follow any 
simple scaling law ( e.g., ~ j 1 - 9 / 2 ) ? corroborating the experimen- 
tal findings. However, both the partial and the sum TCS obey 
different power laws for different incident energy ranges. 

Finally, the present results might be important for the 
future detailed H experiments. 
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